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Abstract.  We  investigate  the  eigenvalue  spacing  distributions  for  ran¬ 
domly  generated  4- regular  Cayley  graphs  on  51/2  (Fp)  by  numerically  cal¬ 
culating  their  spectra.  We  present  strong  evidence  that  the  distributions 
are  Poisson  and  hence  do  not  follow  the  Gaussian  orthogonal  ensemble. 
Among  the  Cayley  graphs  of  SL^iWp)  we  consider  are  the  new  expander 
graphs  recently  discovered  by  Y.  Shalom.  In  additionPwe  use  a  Markov 
chain  method  to  generate  random  4-regular  graphsP  and  observe  that  the 
average  eigenvalue  spacings  are  closely  approximated  by  the  Wigner  sur¬ 


mise. 
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1.  Introduction.  One  of  the  most  remarkable  numerical  discoveries 
of  the  recent  past  is  Odlyzko’s  finding  that  the  spacings  of  the  zeros  of 
the  Riemann  zeta  function  closely  follow  the  Gaussian  unitary  ensemble  of 
random  matrix  theory  [15].  As  a  result  of  this  workFattention  has  turned 
to  the  spacing  distributions  for  the  spectrum  of  other  natural  classes  of 
operatorsF  in  the  hope  of  making  similar  connections  with  other  number 
theoretic  objects. 

One  direction  of  related  work  is  towards  the  analysis  of  the  eigenvalue 
spacing  distribution  for  the  Laplacian  on  different  manifolds;  see  [17]  for  a 
survey  of  many  recent  results  of  this  type  and  an  extensive  bibliography. 
The  motivation  for  our  work  comes  from  the  particular  case  of  interest 
of  the  Laplacian  on  51-2  (Z)  VET  where  IH  denotes  the  hyperbolic  upper 
half  plane.  This  is  the  main  example  of  an  arithmetic  surfaceT  which  is 
a  hyperbolic  surface  given  as  the  quotient  of  the  upper  half  plane  by  an 
arithmetic  subgroup  of  5L2(1R)- 

In  generalFthe  statistical  behavior  of  eigenvalue  spacings  for  natural 
families  of  operators  falls  into  two  main  claissesF  Poisson  and  the  Gaus¬ 
sian  ensemble.  The  Gaussian  orthogonal  ensemble  (GOE)  governs  random 
symmetric  matrices  while  the  Gaussian  unitary  ensemble  (GUE)  applies 
to  random  complex  Hermitian  matrices.  The  density  of  the  spacings  in 
the  Poisson  caseFwhere  the  spacings  are  normalized  to  have  unit  meanFis 
exponentialFe“''Fand  for  the  Gaussian  orthogonal  ensemble  the  density  is 
well-approximated  by  the  Wigner  surmise  [13].  More  generallyF 

Katz  and  Sarnak  [7]  have  recently  investigated  the  eigenvalue  spacings  for 
the  classical  groupsFas  well  as  connections  to  zeta  functions  for  curves  over 
finite  fields.  Diaconis  and  Shahshahani  [3]  have  analyzed  the  eigenvalue  dis¬ 
tribution  for  a  randomly  chosen  matrix  from  classical  groups.  The  physical 
interpretation  of  eigenvalues  is  as  energy  levelsFand  Poisson  behavior  of  the 
spacings  is  usually  thought  of  as  characteristic  of  integrable  systemsFwhile 
GOE  corresponds  to  chaotic  systems  (cf.  [IRIIS]  and  the  many  references 
therein). 

Computations  of  Schmit  [19]  indicate  that  the  spacing  distribution  for 
SL2{'Zi)\M  should  be  Poisson.  More  generallyF  recent  work  of  Luo  and 
Sarnak  [12]  and  Rudnick  and  Sarnak  [16]  indicates  that  this  is  the  case  for 
any  arithmetic  surface.  The  exact  behavior  of  the  spacings  in  this  case  is 
still  an  open  question. 

In  this  paper  we  investigate  the  eigenvalue  spacings  for  the  adjacency 
matrices  of  various  Cayley  graphs  on  5'L2(Fp)Fnaturally  thought  of  as  dis¬ 
crete  approximations  to  the  spectral  behavior  in  the  continuous  setting  of 
S'L2(^)\H.  This  sort  of  analogy  is  suggested  by  the  machinery  developed 
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in  the  successful  application  of  Selberg’s  Theorem  [20]  to  the  discovery  of 
expander  graphs  built  as  Cayley  graphs  of  5^2  (Fp)-  Here  the  main  tool 
is  the  transfer  of  the  lower  bound  on  the  first  positive  eigenvalue  of  the 
Laplacian  to  a  uniform  bound  on  the  first  nonzero  eigenvalue  for  a  family 
of  graphs  obtained  as  qnotients  of  Cayley  graphs  on  SL-^{^).  See  [10]  for 
an  excellent  treatment  of  this  construction  and  a  thorough  bibliography. 

The  results  in  this  paper  are  of  two  types.  On  the  one  hand  we  com¬ 
pute  the  spacing  distributions  for  particular  Cayley  graphs  on  £'1/2  (Fp). 
These  computations  include  generator  sets  to  which  Selberg’s  Theorem 
would  applyFas  well  as  some  sets  to  which  it  would  not  apply.  We  also 
compute  the  distributions  for  random  4-regular  Cayley  graphs  on  S'I/2(Fp) 
(for  a  fixed  prime  p).  In  all  cases  we  find  close  agreement  with  Poisson 
behavior.  For  comparison  we  have  also  computed  the  spectra  for  randomly 
generated  4-regular  graphsF  using  Markov  chain  methods  to  generate  the 
graphs.  These  all  closely  follow  the  GOEPand  this  is  in  agreement  with  the 
extensive  computations  performed  by  JakobsonFMillerFRivin  and  Rudnick 
[5].  Our  computations  on  £'Z/2(Fp)  are  only  made  possible  by  the  use  of 
representation  theory  for  the  groups  51/2  (Fp)  (see  Section  2).  In  briefTwe 
compute  the  spectrum  as  the  union  of  spectra  of  individual  Fourier  trans¬ 
forms  of  the  characteristic  function  of  the  generating  set.  The  second  type 
of  computation  we  carry  out  investigates  the  spectra  of  some  of  these  in¬ 
dividual  transforms.  This  includes  the  spectrum  of  the  expander  graphs 
built  as  the  action  of  51/2 (Fp)  on  the  projective  line.  AgainFin  all  cases  the 
distributions  appear  to  exhibit  Poisson  behavior. 

After  outlining  our  use  of  Fourier  analysis  in  Section  2Fwe  present  our 
computations  in  Section  3.  We  summarize  our  results  as  two  conjectures 
in  Section  4. 

Acknowledgements.  Thanks  to  Peter  Sarnak  for  suggesting  some  of 
these  computationsFand  to  Persi  DiaconisFDmitry  JakobsonFRavi  Kannan 
and  Danny  Sleator  for  helpful  discussions. 

2.  Cayley  graphs  and  Fourier  analysis.  As  in  [819] F  to  analyze 
the  spectrum  of  Cayley  graphs  we  exploit  the  fact  that  the  adjacency  ma¬ 
trix  can  be  viewed  as  the  Fourier  transform  of  the  generators  at  the  right 
regular  representation.  Any  representation  is  equivalent  to  a  direct  sum  of 
irreducible  representations.  Since  we  are  able  to  compute  Fourier  trans¬ 
forms  at  any  irreducible  representation  for  5L2(Fp)Fwe  are  able  to  recover 
the  complete  spectrum  by  only  computing  the  spectrum  of  the  individual 
Fourier  transforms.  For  5L2(Fp)rall  multidimensional  representations  are 
roughly  of  degree  pFversus  (9(p^)Fthe  size  of  the  regular  representation. 

To  say  this  more  preciselyF  let  G  be  a  finite  group  and  let  5  C  G 
generate  G.  The  Cayley  graph  X  —  X{G,S)  for  G  with  respect  to  5  is 
the  undirected  graph  with  vertex  set  equal  to  GF  such  that  there  is  an 
edge  between  a  and  6  in  X  if  and  only  if  as  =  6  for  some  s  G  5  U  5~^. 
EquivalentlyFthe  adjacency  matrix  of  X{G,  S)  has  a  one  in  the  (a,  b)  entry 
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if  and  only  if  as  =  6  for  some  s  G  Let  denote  the  right  regular 

representation  of  GFcomputed  with  respect  to  the  basis  of  delta  functions 
on  G.  Then  it  is  not  difficult  to  see  we  have  the  following  expression  for 
the  adjacency  matrix  of  X(G, .SjFdenoted  F(G,  iS): 

(2.1)  r{G,S)=  Pregi^)- 

sesus-^ 

The  righthand  side  of  (2.1)  is  also  the  Fourier  transform  of  the  character¬ 
istic  function  for  <S  U  5“^. 

Direct  computation  of  the  spectrum  of  the  \G\  x  |G|  matrix  F(G,  S)  re¬ 
quires  0(|Gp)  operations  (cf.  [24]).  For  5L2(Fp)Fsince  |5'L2(Fp)|  =  0{p^)T 
this  means  0(p®)  operations.  This  cost  quickly  becomes  prohibitive.  How- 
everFby  using  the  tools  of  representation  theoryFwe  may  instead  compute 
the  elements  of  an  equivalent  block  diagonal  matrix  and  realize  the  entire 
spectrum  as  the  union  of  the  spectra  of  the  blocks. 

More  preciselyFrepresentation  theory  gives  a  simultaneous  block  diag- 
onalization  of  the  matrices  p^egi^) 
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0 
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where  pi, . . . ,  ph  are  a  complete  set  of  irreducible  matrix  representations  of 
G  with  deg(/9i)  copies  of  pi{s)  on  the  diagonal.  ConsequentlyF 

(2.4)  spectrum[X {G ,  S))  =  [J  spectrum  I  Pi{s)  1  ■ 

»  =  1  VsgSUS-i  / 

The  degree  of  the  largest  irreducible  representation  of  G  is  bounded  above 
by  |G|^/^  (cf.  [21]).  ThusFby  using  the  basis  which  gives  (2.2)  and  (2.3)Fwe 
are  able  to  reduce  the  computation  from  0(|G|^)  operations  to  a  potentially 
more  manageable  0(|G|  •  (maxj  deg(/7j)))  <  OdG^'^^)  operations. 

For  G  =  5L2(Fp)Fthe  irreducible  representations  of  5'L2(Fp)  occur  in 
two  familiesFthe  discrete  series  and  principal  series.  The  distinction  de¬ 
pends  upon  the  restriction  of  an  irreducible  representation  to  the  Borel 
subgroup  B  <  SL2{Wp)  of  upper  triangular  matrices.  An  irreducible  repre¬ 
sentation  of  51/2 (Fp)  is  said  to  be  from  the  principal  series  if  its  restriction 
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to  B  contains  the  trivial  representation.  OtherwiseFit  is  said  to  be  from 
the  discrete  series.  The  principal  series  representations  occur  as  compo¬ 
nents  of  induced  1-dimensional  representations  from  SF which  are  almost 
always  irreducible.  This  gives  the  trivial  representationFone  representation 
of  degree  pFtwo  representations  of  degree  {p+  l)/2Fand  [p  —  3)/2  repre¬ 
sentations  of  degree  p+1.  The  discrete  series  is  less  easily  explainedFbut 
suffice  it  to  say  here  that  the  representations  are  in  close  correspondence 
to  the  characters  of  the  non-split  torus  in  SL2{Wp)  (cf.  [14]FCh.  2FSection 
5).  There  are  two  such  representations  of  degree  [p  —  l)/2  and  (p  —  l)/2 
representations  of  degree  p  —  1. 

Explicit  representations  are  needed  to  apply  (2.2)  and  (2.3).  These 
can  be  found  in  [14]  and  [23]  and  are  the  basis  of  our  implementation 
(cf.  [8]).  Knowledge  of  the  representations  of  51/2  (Fp)  gives  the  irreducible 
representations  of  PSL2{¥p)-  More  preciselyFif  p  is  an  irreducible  matrix 
representation  of  51/2  (Fp)  and  —I  is  in  the  kernel  of  pF  where  I  —  (JJ)F 
then  p  is  constant  on  cosets  5Z/2(Fp)/{±f}  and  as  such  gives  an  irreducible 
representation  of  P5Z/2(Fp)-  Under  this  identificationFthe  set  {p  \  {±J}  C 
ker(p)}  gives  a  complete  set  of  inequivalent  irreducible  representations  of 
PSL2{¥p).  We  use  this  correspondence  when  calculating  the  spectra  for 
the  Ramanujan  graphs  constructed  by  LubotzkyF Phillips  and  Sarnak  in 
Section  3. 

Remark.  In  the  case  in  which  <S  U  is  a  union  of  conjugacy  of 
classes  of  GFthe  adjacency  matrix  can  be  completely  diagonalized  and  the 
eigenvalues  can  be  computed  as  certain  character  sums  over  G  [2]. 

3.  Numerical  evidence.  The  main  results  of  this  note  are  computa¬ 
tions  of  the  spacings  for  various  Cayley  graphs  on  5L2(Fp)Fmade  possible 
by  the  techniques  outlined  in  Section  2.  As  we  stated  in  the  introductionF 
our  computations  are  of  two  types: 

(1)  The  computation  of  spacing  distributions  for  particular  2-generator 
Cayley  graphs  on  5L2(Fp). 

(2)  The  computation  of  spacing  distributions  for  particular  Fourier 
transforms  for  2- generator  Cayley  graphs  on  SL2{Wp)- 

Within  these  computations  we  distinguish  among  various  kinds  of  gen¬ 
erators;  global  or  Selberg  generators  are  generators  for  families  of  Cayley 
graphs  obtained  as  projections  of  a  single  Cayley  graph  on  SL2(’^)',  non- 
Selberg  generators  are  an  infinite  family  of  generators  for  5Z/2(Fp)  as  p  — >  oo 
that  are  not  the  projection  of  a  single  set  of  generators  for  5^2  (^);  fan¬ 
dom  generators  are  a  pair  of  generators  for  SL2{¥p)  generated  by  a  simple 
randomized  algorithm.  In  general Frandom  generators  are  non-Selberg. 

For  a  generating  subset  S  C  SL2{Fp)  we  denote  the  eigenvalues  (with¬ 
out  multiplicities)  of  F(G,>S)  as  Aq  >  Ai  >  •  •  •  >  AjvFand  let  P{S)  denote 
the  “empirical”  cumulative  distribution  function  for  the  eigenvalue  spac- 
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ingsFso  that 


(3.1)  PiS)  =  ^f^[\j^r-Xj<S] 

i=i 

where  [a  <  6]  is  one  if  a  <  6  is  true  and  is  zero  otherwise.  We  assume  the 
eigenvalues  are  normalized  so  that  the  spacings  have  mean  one: 

1  ^ 

(3.2)  = 

j=i 

Besides  the  computations  of  the  various  P(S)rwe  also  include  some 
new  data  on  a  particularly  interesting  non-Selberg  generating  pair  recently 
discovered  by  Y.  Shalom  [22].  These  turn  out  to  be  (in  terms  of  the  nu¬ 
merical  analysis  of  the  second-largest  eigenvalue)  among  the  best  expanders 
built  as  Cayley  graphs  of  5Z/2(Fp)  discovered  to  date. 

3.1.  Random  graphs.  In  this  section  we  present  sample  results  from 
our  calculations  of  the  spacing  distributions  of  random  4-regular  graphsF 
where  the  curves  show  GOE  behavior  (see  also  [5])Fand  4-regular  Cayley 
graphs  of  5Z/2(Fp)  for  randomly  chosen  generating  pairs Fwhere  the  behavior 
is  Poisson. 

To  choose  generating  pairs  of  5L2(Fp)  uniformly  at  randomFwe  use  the 
algorithm  described  in  [8] .  This  algorithm  first  selects  two  group  elements 
0,6  €  SL2{]Fp)  uniformly  at  random  using  the  Bruhat  decomposition  of 
>S'L2(Fp)Fand  then  checks  whether  {a,  6}  generates  the  group  by  verifying 
that  {7r(a),  7r(6)}  does  not  generate  one  of  the  six  possible  subgroups  of 
PiS'L2(Fp) Fwhere  tt  :  5L2(Fp)  — t  PSL2{Fp)  is  the  natural  projection.  We 
refer  to  [8]  for  details. 

We  generated  random  Cayley  graphs  using  this  algorithm  for  p  >  150F 
and  observed  that  all  of  the  graphs  closely  followed  the  exponential  distri¬ 
bution.  In  Figure  3.1  we  show  the  cumulative  distribution  function  for  a 
typical  example. 

Remark.  In  Figure  3.1Fcis  in  all  of  the  cumulative  distributions  that 
we  presentFwe  have  omitted  the  spectral  gap  k  —  \i  from  the  calculations. 
AsymptoticallyFas  p  ooFthis  gap  does  not  contribute  to  the  cdfFalthough 
it  changes  the  mean  of  the  distribution. 

To  generate  random  4-regular  graphs  (not  Cayley  graphs)  we  used  a 
Markov  chain  method.  In  generalFthe  states  of  the  Markov  chain  are  the 
(labeled)  A:-regular  graphsF  and  two  graphs  can  be  connected  by  a  single 
step  of  the  random  walk  if  and  only  if  the  symmetric  difference  of  their 
edges  is  a  cycle  of  length  4.  The  random  walk  can  be  described  in  terms  of 
the  incidence  matrices  of  the  graphs.  Recall  that  the  incidence  matrix  of  a 
graph  F(£',  V)  is  the  jlAl  x  \E\  matrix  /(F)  where  the  column  corresponding 
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Fig.  3.1,  The  cumulative  distribution  function  P(S)  for  the  eigenvalue  spacings  of  the 
single,  randomly  chosen  generating  pair  a  =  (  Isa)’  **  =  (llo  124) 

The  curve  does  not  include  the  spacing  between  the  first  and  second  eigenvalues,  which 
in  this  case  was  1  —  0.879090.  The  dashed  line  is  the  curve  1  —  e~^ . 


to  edge  (i,  j)  has  a  1  in  the  and  j*'*  rowFand  O’s  elsewhere.  If  I  is  the  state 
of  the  random  walkFtwo  rows  1  <  i,  j  <  |  V|  and  two  columns  1  <  k,l  <  \E\ 
are  chosen  uniformly  at  random.  If  lik  =  Ijt  =  1  and  /,•;  =  Ijk  =  0  then 
the  chain  moves  to  the  state  with  7,^  =  Iji  —  0  and  In  =  Ijk  =  \  unless 
a  double  edge  is  formed  by  doing  so.  SimilarlyFif  lik  =  Iji  =  0  and 
In  —  Iji,  =  1  then  the  chain  moves  to  the  state  with  lik  =  i^ji  =  ^  find 
In  =  Ijj,  =  OFagain  unless  a  double  edge  would  be  formed  by  this  move. 
In  all  other  cases  the  walk  remains  in  the  same  state. 

It  is  proved  in  [6]  that  this  random  walk  is  rapidly  mixingFusing  the 
technique  of  canonical  paths  to  estimate  conductance.  While  this  results 
in  a  polynomial  algorithm  in  the  size  of  the  graphsFthe  exponent  is  too 
large  to  enable  this  result  to  yield  a  stopping  criterion  for  the  graphs  we 
generate.  As  a  matter  of  practicalityFwe  simply  run  the  chain  for  a  large 
(10®)  number  of  steps  to  generate  each  graph.  After  stopping  the  chainF 
we  test  to  make  sure  the  graph  is  connected  (it  is  with  high  probability) . 

Since  the  Cayley  graphs  for  57/2 (F157)  yield  on  the  order  of  25ID00 
eigenvalues  and  19K)00  intervalsFwe  generated  10  random  4-regular  graphs 
on  2ID00  vertices  and  averaged  the  intervals  to  obtain  comparable  statistics. 
The  resulting  cumulative  distribution  function  is  shown  in  Figure  3.2Fwhere 
it  is  compared  with  the  Wigner  surmise. 
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Fig.  3.2.  The  cumulative  distribution  function  P{S)  for  the  eigenvalue  spacings  of  the 
10  4-regular  graphs  on  2,600  vertices,  generated  by  running  the  Markov  chain  for  10® 
steps.  The  average  value  of  Xi  was  0.86SS89.  The  dashed  line  is  the  cdf  for  the  Wigner 
surmise,  1  —  exp(— 


3.2.  Explicit  generators.  When  we  computed  the  spacing  distribu¬ 
tions  for  both  Selberg  and  non-Selberg  type  generatorsF and  for  the  gener¬ 
ators  recently  discovered  by  Shalom  [22]rwe  found  that  they  followed  the 
Poisson  behavior  very  closely.  The  generators  for  5^2  (Fp)  we  used  were  the 
following: 


(3.3)  {(J  ;)}  Selherg 

(3.4)  IQ  l)’(3  l)}  non-Selberg 

(3-5)  iTui)} 


where  w  is  a  primitive  cube  root  of  unity  {mod  p).  For  p  =  199rthe 
distribution  for  Shalom’s  generators  is  shown  in  Figure  3.3.  The  curves  for 
the  other  generators  are  very  similar. 

The  second-largest  eigenvalue  Ai  is  shown  for  a  few  of  the  small  primes 
p  =  1  {mod  4),  in  Table  3,2.  These  numbers  indicate  thatFexcept  for  the 
known  Ramanujan  graphsF these  generators  are  perhaps  the  best  explicit 
4-regular  graphs  for  57/2 (Fp)  that  have  been  obtained.  For  P57/2(Fp)Fthe 
explicit  Ramanujan  graphs  of  LubotzkyF Phillips  and  Sarnak  have  better 
separation.  For  comparisonFthe  second-largest  eigenvalue  for  the  LPS  gen- 
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Fig.  3.3.  The  cumulative  distribution  Junction  P(S)  for  generators  (S.S)  with  p  =  199, 
so  a  =  ^i*)>  ^  —  ( 107  uir)'  second-largest  eigenvalue  is  Ai  «  0.886048. 


erators 

(3.6) 


1  /  i  1  ±  i  \  _1_  f  i  -1  ±  A 

^  V  “1  ±  *  -*  /  ’  V 1  ±  / 


where  i  =  t/— IFare  shown  for  several  values  of/?  =  3  [mod  4)  in  Table  3.2. 

Not  all  of  the  spacing  distributions  that  we  observed  for  5L2(Fp)  were 
so  closely  Poisson.  As  an  examplePthe  plot  in  Figure  3.4  shows  the  spacing 
distribution  for  the  Selberg  type  generators 


for  p  =  199.  We  believe  that  the  discrepancy  can  be  attributed  to  the 
expansion  properties  of  the  graph.  (The  difference  cannot  be  attributed 
to  gaps  in  the  spectrumF  since  other  than  the  gap  4  —  AiFthe  spectrum 
is  “continuous.”)  The  second-largest  eigenvalue  for  this  graph  is  Ai  » 
0.977554. 

3.3.  Individual  transforms.  The  data  presented  above  strongly  in¬ 
dicates  that  the  average  spacing  associated  with  the  Fourier  transforms  of 
the  delta  function  supported  on  the  generating  set  is  asymptotically  Pois¬ 
son.  In  this  section  we  investigate  the  behavior  of  the  individual  transformsF 
and  present  sample  calculations  for  the  LPS  graphs  [11E18].  Of  partic¬ 
ular  interest  is  the  transform  at  the  principal  series  representation  induced 
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p 

Ai 

P 

Ai 

13 

0.921721 

109 

0.908280 

19 

0.864292 

127 

0.886755 

31 

0.879056 

139 

0.887515 

37 

0.894603 

151 

0.880739 

43 

0.894659 

157 

0.888702 

61 

0.882493 

163 

0.881105 

67 

0.890386 

181 

0.890412 

73 

0.903236 

193 

0.887729 

79 

0.881429 

199 

0.886048 

97 

0.899836 

211 

0.881852 

103 

0.883868 

223 

0.882423 

Table  3.1 

Ai  for  a  sequence  of  Shalom’s  expanders. 


P 

Ai 

13 

0.832880 

37 

0.863086 

61 

0.865375 

73 

0.862093 

97 

0.864023 

109 

0.863180 

157 

0.861790 

181 

0.863598 

193 

0.862532 

229 

0.865491 

241 

0.864479 

Table  3.2 

Ai  for  4-regular  Ramanujan  graphs  (3.6)  on  PSL2{Wp).  The  Ramanujan  bound  is 
V3/2  w  0.8660254. 
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Fig,  3.4.  P{S)  with  generators  a=(j  })>^  ~  (-i  i)  ‘^■^2(Fi99)-  The  second- 
largest  eigenvalue  is  Ai  fts  0.977554. 

from  the  identityF since  this  is  associated  with  a  graph  on  the  projective 
line. 

To  explainFwe  recall  that  SL2{Wq)  acts  on  the  projective  line  P^(F9)  = 
{0, 1, . . . ,  9  —  1,  oo}  by  fractional  linear  transformations: 

I'  a  b\  _  auj  +  b 

\  c  d  J  coj  d' 

Let  B  denote  the  Borel  subgroup  |  ^  “  JLi  j  j-  and  let  the  matrices 

Su  =  ^  ^  ,  M  =  0, 1, . . . ,  g  -  1  and  Sqo  = 

be  a  fixed  set  of  coset  representatives  of  B\S'L2(F5)Fwhich  we  identify  with 
F^(Fq)-  In  [8]  we  constructed  the  principal  series  representations  of  5^2  (Fj) 
by  inducing  characters  from  BT  and  expressing  these  matrices  in  terms  of 
their  action  on  {«„}.  Very  brieflyFif  V'  :  F,  — t  C  is  a  character  of  BFthen 
the  induced  representation  =  ibt B  is  given  as 

(3-8)  P,pig)su  =  i’ifiu,g))sg.s, 

for  some  function  f{u,s)  G  F^-  In  particularF (J5 (pi )  is  the  adjacency 
matrix  F(P^(Fq),  <5);  we  refer  to  [8]  for  details.  If  T{PSL2{Wq),S)  is  a  good 
expanderFthen  F(P^(F5),  >5)  is  also  a  good  expander  since  its  spectral  gap  is 
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Fig.  3.5.  P{S)  for  individual  Fourier  transforms,  for  LPS  generators  with  p  =  37  and 
q  =  1009.  Top  row:  two  principal  series  representations,  action  on  projective  line  on 
the  left.  Bottom  row:  two  discrete  series  representations. 

at  least  as  large.  In  particularPthe  explicit  Ramanujan  graphs  for  PSL2{Wq) 
restrict  to  give  Ramanujan  graphs  on  P^(Fg)  [18]. 

From  a  computational  point  of  viewFthe  discrete  series  representations 
are  more  expensive  to  construct.  HoweverFour  experience  has  shown  that 
the  spectral  properties  of  the  principal  and  discrete  series  representations 
are  almost  identical.  In  Figure  3.5  the  spacing  distribution  is  shown  for 
the  38-regular  LPS  graphs  P5L2(Fioo9)  two  principal  series 

and  two  discrete  series  representations. 

4.  Summary.  We  have  computed  the  cdf  for  the  eigenvalue  spacings 
of  a  variety  of  Cayley  graphs  on  SL2{Fp)-  When  the  graphs  are  good  ex- 
pandersFthe  distributions  are  PoissonFeven  on  the  level  of  the  individual 
Fourier  transforms.  The  data  suggests  and  supports  the  following  conjec¬ 
ture. 

Conjecture  1.  Let  Xp  —  X{PSL2{¥p),S)  be  a  family  of  k-regular 
Cayley  graphs  such  that  Ai(rj,)  — >  2y/k  —  1/k  as  p  cx>.  Then  asymptoti¬ 
cally  as  p  —>  oo,  the  distributions  of  the  eigenvalue  spacings  for  Xp  and  its 
individual  Fourier  transforms  Ss{p)  are  Poisson. 

AdditionallyF  our  randomly  generated  2-generator  Cayley  graphs  on 
S'L2(Fp)  also  exhibit  Poisson  behavior.  In  light  of  our  previous  experiments 
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[819]  which  suggest  that  a  randomly  chosen  2-generator  Cayley  graph  on 
SL2{Wp)  is  a  good  expanderF this  too  suggests  a  relationship  between  ex¬ 
pansion  and  Poisson  spacing.  In  this  spirit  we  anticipate  a  “central  limit 
theorem”  for  random  2-generator  Cayley  graphs  on  5L2(Fp)- 

Conjecture  2.  For  e  >  0,  as  p  oo  the  probability  of  choosing  a  pair 
of  generators  S  for  SL2{Wp)  with  ||P(5')  —  (1  —  e“®)||  >  £  goes  to  zero. 

Availability  of  Software.  All  of  our  computations  were  performed 
on  HP  735/125  and  DEC  3000  Model  600  Alpha  workstationsFusing  soft¬ 
ware  written  in  the  C  language.  The  code  will  be  available  via  the  web 
page  WWW .  cs  .  dartmouth .  edu/~rockmore/GFT. 
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